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Abstract 

The competition among spin glass (SG), antiferromagnetism (AF) 
and local pairing superconductivity (PAIR) is studied in a two-sublattice 
fermionic Ising spin glass model with a local BCS pairing interaction in 
the presence of an applied magnetic transverse field F. In the present 
approach, spins in difi'erent sublattices interact with a Gaussian random 
coupling with an antiferromagnetic mean Jo and standard deviation J. 
The problem is formulated in the path integral formalism in which spin 
operators are represented by bilinear combinations of Grassmann vari- 
ables. The saddle-point Grand Canonical potential is obtained within the 
static approximation and the replica symmetric ansatz. The results are 
analysed in phase diagrams in which the AF and the SG phases can occur 
for small g {g is the strength of the local superconductor coupling written 
in units of J), while the PAIR phase appears as unique solution for large 
g. However, there is a complex line transition separating the PAIR phase 
from the others. It is second order at high temperature that ends in a 
tricritical point. The quantum fluctuations affect deeply the transition 
lines and the tricritical point due to the presence of F. 



1 Introduction 

It is now well-established that strongly correlated systems such as heavy fermions 
(HF) P and high temperature superconductors (HTSC) [2], upon doping, can 
present magnetic order or superconductivity. The complexity involved in such 
physical systems as, for instance, the existence of Non-Fermi Liquid (NFL) be- 
haviour, has given rise to new theoretical approaches. In particular, recent works 
have proposed that the presence of disorder can aflFect these strongly correlated 
systems (see j3j and references therein) being even source of NFL behaviour, 
for example, in HF The presence of disorder can also induce frustration 

which, in fact, has been found in several HF 00111111 and HTSC HOIIIII phys- 
ical systems. Some theories have been investigating whether or not spin glass 
(SG) phase can be found in models designed to study certain aspects of HF or 
HTSC systems ^lElEl- For instance, the existence of a SG solution has 
been demonstrated in the Kondo lattice model |15[ll6j . Nevertheless, relatively 
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little consideration 17 has been given in order to obtain the behaviour of the 
transition temperatures for these strongly correlated physical systems and, thus, 
to mimic the phase boundaries of their global phase diagram which includes an- 
tiferromagnetism (AF), SG, superconductivity and the possible presence of a 
Quantum Critical Point (QCP). 

Regarding the experimental scenario, there are some similarities in the global 
phase diagrams between some HF and HTSC, although the microscopical mech- 
anisms involved in such systems are different |18| . For instance, the HTSC 
compound Yi-xCax:Ba2Cu30e has a phase diagram temperature T versus the 
hole concentration psh [101 which displays an AF ordering for low psh- The 
respective Neel temperature Tn decreases for < Psh < 0.035 until the onset 
of a second transition at Tf in which is found the superposition of frustration 
with a preformed AF background. For psh > 0.035, there is another transition 
from the previous mixed state to a pure SG. The SG transition temperature 
Tg also decreases with the increase of Psh until the onset of superconductivity 
(SC). However, the SG regime is still found into the SC region. For instance, for 
0.06 < Psh < 0.10 there is still traces of spin freezing. Finally, for large values of 
Psh, the superconductivity is dominant. From the side of the HF systems, the 
example is the compound Ui-xLaxPd2Al3. When the doping of La is increased, 
the corresponding phase diagram shows the AF ordering to be replaced by a SG 
state in the region 0.25 < a; < 0.65. At a; = 1, the system is a superconductor. 
The Neel temperature Tn is sharply decreased from T = 14.6K until T — 2.6K 
when X increases, while the subsequent SG transition temperature Tg drops to 
a QCP at a; ~ 0.8. In the intermediated doping region between the QCP and 
X = 1, a, NFL behaviour is observed. 

In the last years, several works have been studying the competition between 
SG and pairing formation in a formulation where the spins are represented as bi- 
linear combination of fermionic creation and destruction operators. |19II2(J[|^ . 
The model used in such approaches is composed of a random Gaussian cou- 
pling between localized spins together with pairing interaction in the real space. 
In fact, it can be shown that both terms of this simple model have the same 
origin. They can be derived by eliminating the conduction electrons, in sec- 
ond order of perturbation (see Appendix in Ref. 1 1 9J ) ■ from an earlier model 
introduced to treat conventional superconductor doped with magnetic impuri- 
ties 122]. A saddle-point solution for the respective grand-canonical potential 
has been obtained within functional integral formalism for the Ising |19j and 
the Heisenberg version of the model in the half-filled situation using the 
static approximation |23j and the replica symmetry ansatz. For the Ising case, 
the phase diagram temperature versus the strength of the pairing coupling g 
(in units of J that is the variance of the Gaussian random distributed spin-spin 
coupling ) displays two phase boundaries. For lower it is found a second order 
line transition between paramagnetism (NP) and SG at Tg = 0.95J. For large 
g, there is a complex line transition Ti(g) separating the PAIR phase (where 
pair formation is found) from the SG and NP phases. The line transition Ti{g) 
is a second order for high temperature when the boundary is between NP and 
the PAIR phase. However, it becomes first order at lower temperature. The 
corresponding tricritical point Ttrict has been located on Ti (g) above of Tg . As 
consequence, the boundary between SG and the PAIR phase is entirely first or- 
der. Weak hopping corrections done elsewhere j2Il demonstrated that the phase 
boundaries described previously with no hopping are essentially maintained. 
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The model in Ref. has two important sliortcomings. The first one is the 
lack of quantum spin flipping mechanism which would be able to suppress the 
transition temperatures leading them to a QCP. Even the Heisenberg extension 
of the problem [21] has been unable to produce a QCP. This particular weakness 
has been corrected in Ref. by the addition of a transverse field F in the 
Ising version of the model. The presence of F has changed the behaviour of both 
transition temperatures Tg and As long F increases, the first one moves 

downwards in the direction of a QCP while the second one is displaced. There- 
fore, F has also suppressed the PAIR solution in the sense that it is necessary 
larger values of the pairing coupling strength g to find the PAIR solution |24j . 
The tricritical point Ttrict is affected by the presence of F. The transverse field 
moves up Ttrict which enlarges the first order portion of Ti{g). Finally, it has 
been proposed a relationship between F and g (J is kept constant) based on the 
argument that the pairing and RKKY interaction have the same origin in the 
derivation of the model [l^j. As consequence, the effects described previously 
are superposed in a single phase diagram T versus g due to the increase of g, 
hence F. It shows Tg decreasing towards a QCP at g = gc, then a PAIR phase 
can be found a,t g > gc with Ttrict located at higher values of T and g than the 
case F = 0. 

Nevertheless, the model used in Refs. EH |^ has a second shortcom- 
ing, it is unable to produce an AF solution. Thus, the model is useless if one 
is trying to study the phase boundaries which include SG, PAIR phase and 
also AF. However, quite recently the competition between AF and SG has been 
analysed in a disordered two-sublattice fermionic spin model. The model is a 
Gaussian random coupling with an antiferromagnetic mean Jo and standard 
deviation J between spins in distinct sublattices with the presence of a trans- 
verse F and parallel H magnetic fields .25. . In fact, it is the fermionic version 
of the model introduced by Korenblit and Shender (KS) [201 which is used to 
study the competition between AF and SG with classical Ising variables. This 
classical two-sublattice model has itself unexpected effects as compared with the 
classical single lattice SG Ising model |2ZI. For example, opposite solutions are 
enforced by the degree of frustration {Jo)~'^ and H (given in units of J). When 
degree of frustration decreases the AF solution is favored, while the field H can 
eventually enhance the frustration in a certain range f^. This last effect is 
related with the asymmetry between the two-sublattice due the coupling with 
H. The presence of the F in the fermionic version of KS model has introduced 
important differences if compared with its classical counterparts as long there is 
a competing mechanism associated with Jq, H and F |25'. For instance, F sup- 
presses the magnetic orders leading their critical temperatures to QCPs, while 
H can favour frustration at the same time that it destroys the AF phase. 

Therefore, the purpose of the present work is to study the competition among 
AF, SG and the PAIR phase using the fermionic Ising KS model with a local 
pairing interaction in each sublaticce in the presence of a transverse field F. 
Particularly, the focus is to describe the behaviour of the possible transition 
temperatures present in the problem. We follow the same approach used in 
Refs. ^Jl^nil^. The partition function is obtained in the functional inte- 
gral formalism where the spin operators are given as bilinear combinations of 
Grassmann fields. The static approximation (SA) and the replica symmetry 
(RS) ansatz are used to calculate the saddle-point Grand canonical potential. 
Particularly, we extend for the present two-sublattice problem a procedure that 
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mixes Nambu matrices and spinors as already introduced in Ref. li.. The 
stability of the RS solution is also investigated. Surely, the simple model used 
in the present work is not suited to describe the extremely complicate physics 
present in HF as well as in HTSC systems. However, it can be, at least, useful 
to mimic general features of a phase diagram in which AF, SG, pairing coupling 
and a quantum spin flipping mechanism are present. 

The use of SA and RS ansatz deserves some remarks. It is well known the 
SA is not adequate to describe the low temperature behaviour of the spin-spin 
correlation function |2S|- However, the justification for the use of the SA in 
this work is based on the fact that our interest is mainly to obtain the possible 
transition temperatures associated with AF, SG and PAIR competition. The 
analysis of the quantum rotor model in the Af — > cxi indicates that the critical 
line can be obtained from the zero frequency mode 28' . 

Although the focus in this work is to obtain the phase boundaries of the 
fermionic version of KS model with a pairing interaction in the presence of 
a transverse field F, the thermodynamics is not the only source of valuable 
information. In the case of Ising SG fermionic model, other physical quantities 
can have an interesting behaviour. For example, the density of states (DOS) 
and, hence, the local Green's function are affected at low temperature by the 
replica symmetry breaking . In particular, there is a presence of a pseudogap. 
Recently, a mapping between the one-lattice fermionic Ising SG and the classical 
Ghatak-Sherrington model [3T] has demonstrated that the true origin of such 
effects in the DOS are, in fact, classical 122]. However, for the present model 
we can speculate if there is such kind of mapping due to the presence of the 
transverse field. In that sense, this work can be also thought as a first step 
towards the understanding of the problem given by the model introduced here 
from a many-body perspective since the thermodynamics is well understood, at 
least, at mean field level. 

This paper presents the following structure. In section 2, we derive the 
saddle-point Grand Canonical potential and the set of equations for the order 
parameters which is enlarged when compared with a single lattice Ising model 
studied in Ref. |24| . In section 3, we solve the order parameter equations. 
In order to capture properly the competition among the phases present in the 
problem, we build up phase diagrams T versus g for several values of F and 
Jo given in units of J where Jq and J are the mean and standard deviation of 
the random Gaussian spin-spin interlattice coupling, respectively. On the other 
hand, the effects of the transverse field F are better shown in a phase diagram 
T versus F for constant values of g and Jo- We also obtain a phase diagram 
T versus g in which F and Jq are related with g based on the same arguments 
proposed in Ref. j24| . This procedure mixes the effects of both parameters in 
the the phase diagram. In the last section, we make our conclusions. 

2 General Formulation 

The model studied here is composed by interlattice Gaussian random spin-spin 
interaction j26| . a intrasite local BCS pairing interaction (which favors double 
occupation of sites in each sublattice) with a transverse magnetic field applied 
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r. Therefore, the hamiltonian is given by: 

N 

^ = -Yl ■^'-3" - 2r ^ ^ 

iajb P jp = l 

"^im^lT^li^Jpi^JpT (1) 
p ^pjp 

where the sums over ip (jp) are run over the N sites of each sublattice p {p = a 
or b). The exchange interaction Ji^j^ is an independent random variable with 
the following Gaussian distribution: 



N 



64J2 
N 



(2) 



The spin operators in Eq. Q are defined in terms of fermion operators: 



S^P = 2 [^»PT - : SI = - 



where nij,a gives the number of fermions at site ip with spin projection a =t or i- 
c| ^ and Ci a are the fermion creation and annihilation operators, respectively. 

The problem is formulated in a path integral formalism in which the spin op- 
erators are represented as anticommuting Grassmann fields (0*, <^). Therefore, 
the Grand canonical partition function is given by: 

Z = y i?[0*0]exp[ A] (4) 

with the action 

^= rrf^{EE[C^M(7r-^)'^»^-(^)] -i?(r(r),0(r))}, (5) 

^0 P,o ^, 

[3 = 1/T (T is the temperature), r is a complex time and /i is the chemical 
potential. The Fourier decomposition of the time-dependent quantities is em- 
ployed in Eq. (|SJ). The action can be write ds, A = At -\- Asg + Abcs with: 

= E E E ('^) + f^^' + /^ra-]^, ('^)' (6) 

^«G = EE^'^^<-j-^^t(^')^j;(-^')' (7) 
ABC5 = f EEE<(-K(-) (9) 

where Pip^u) ) = 5I],„ (/ii^x (-a;)0ip| (cj +uj), a" (v = x, y or z) denotes the Pauh 
matrices, 0j = (0* 1(1^)) is a Grassmann spinor, and lu — (2to + l)7r 

and Lo — m-K {m — 0, ±1, • • • ) are the Matsubara's frequencies. 
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The grand canonical potential is obtained within the static approximation 
which considers uj = in Eqs. 0-® j^UHj. The configurational averaged 
thermodynamic potential per site is obtained with the use of the replica method: 
/3r2 — —5^ Hni((Z")j.^ — l)/n where the repHcated partition function is: 



{n.h, = I X{D{<i^*'^r)eMAf+A^iG+A%cs] (10) 

where A" is given by Eq. JHl with a sub-index a, 

= E[^(Ec^>)' - ^E^^^] (11) 



^Bcs = ^ E E E Ei"^(-k^i"(-)]^ (12) 

a.p v—x^y ip w 

with the replica index a running from 1 to n. In Eq. (|12|l . it is introduced the 
Nambu matrices /"^(w) and ^'"(w) in which ^'"^(w) = ((/)*f (w) (l)f^{-uj)). 

Eq. Hll(l can be rearranged reviewing the sums over different sublattices 
by square sums over the same sublattice. The replicated partition function is 
then linearized by using Hubbard-Stratonovich transformations. It inserts the 
replica-dependent auxiliary fields q^^ , nip, 77^^ and 77" ^ in Eq. ifTUI) . The 
Gaussian integrals over these fields have been exactly performed in the thermo- 
dynamic limit by the steepest descent method. Therefore, the Grand canonical 
potential is: 

Z{n)/N = /3gE Kl' " ^"^o E « + ^'^^ E ^'^-^^f " ^^^"o.p^'o.p (13) 

a.p a af3 



where 77^ = 77^ p -f i-q'^ ^ and 



- 2/3 Jo E ™? + At,p + /35 E "!Z;t"] (14) 

with the matrix 77^ = rj'^^^a^ + 77" ^ct^. The fields q^^ , nip and \r)p\ in Eq. 

Ijll-ill are given by saddle-point equations, in which qp^ is related with the spin 
glass order parameter, nip is the magnetization of the sublattice p, and [77^1 
is an order parameter that indicates long range order where there is double 
occupation of sites in sublattice p. 

In the present work, it is assumed the replica symmetric ansatz, which con- 
siders = qp for all a ^ /?, g^" = = Xp + Qp, rrip = nip, and 77^ = 77p 
for all a. The physical quantity Pxp is the static susceptibility when Jq = 0. 
The sums over replica indices are performed. It produces quadratic terms in 
Eq. (|14|l that are linearized introducing new auxiliary fields in Eq. H13|l . The 
resulting functional Grassmann integral is an exponential that sums quadratic 
forms of spinors and Nambu matrices. In order to perform the integral over the 
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Grassmann fields, it can be used a matrix that mixes elements of spinors and 
Nambu matrices, such as: 



where Dx — dxe^^'^ / ^/2Tr {x ~ or Zp), 



with 



= [0ptH 't>*pi{w) (t>ip{-w) (t>p^{-w)\ 



(15) 

(16) 
(17) 



G-\uj) 







\ 



\ 







-Pg-n*p 



(18) 



-/3r itj-C" / 



and — (3^p ± /3/ip. The internal field hp = J{^2q^,Zp + ^J^Xp'^p) — ^o'Tip' , 
which acts on the sublattice p, depends on the order parameters of sublattice 

The functional integral in Eq. (|16|l and the sum over the Matsubara's fre- 
quencies can be performed: 



I{zp,^p) = coshpJfj,^ + g^rj^ + cosh^^Ap 



(19) 



with Ap = hp + T^. This result and Eq. H13() are used to express the thermo- 
dynamic potential as: 



2f3n = -pjomarrib + f3^J^iXaXb + XaQb + XbQa) 

/OO 
i^Zplni^p(zp)+ln4 
-OO 



p—a,b ^ 



where 



Kp{zp) = i cosh Pgr]p + / cosh /3y^A^ 

\ J — OO 



(20) 



(21) 



with yti = to ensure the half-filling situation. The order parameters are given 
by the extreme condition of Eq. (|20|l : 



Vp 



Dzp 

_ 1 
" 2 



Kp{zp) 
smh{P grip) 



Dzp — - 

-OO ^p\^p) 



(22) 
(23) 
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Figure 1: Phase diagrams as a function oi T/J and pairing coupling g/J for 
Jo = 1.5 J and for two values of F/J: F = and F — IJ. The full lines indicate 
second-order transitions while the dashed lines indicate first-order transitions. 



Xp 



DZr, 



/"oo^Cp^sinh/jyA; 

Y P 

Kp{zp) 



Coo w ^ cosh /j^a; 

Dz„ 



(24) 



(25) 



Kp{zp) 

The stability of replica symmetric solution is analysed by Almeida- Thouless 
eigenvalue Xat'- 



/oo 
Dz 
-oo 



{Kp{zp)y 



(26) 



where 



1 9^ 

Ip{zp) = Kp{zp) I D^p — — cosh/J^A^ 



(/ D^p^ainhp^p) 



(27) 



3 Phase diagrams 

Numerical investigations of the order parameter equations (|22l25ll allow us to 
find three kinds of solutions. The SG solution corresponds to qa = qt 
(with rua ~ TTib = 0, T]a = Vb = 0) whilc the AF solution is nia = —rrib ^ 
(with qa ~ <lb ^ 0, rja — r?6 ~ 0). The spin pairing solution (PAIR phase) 
corresponds to rja and r/;, different from zero while the rest of order parameters 
is zero. The instability of the replica symmetry (RS) solution of the SG is also 
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Figure 2: Phase diagrams as a function of T/J and pairing coupling g/J for 
Jo = 1.7J and for two values of V / J: F = and F = IJ. It is used the same 
convention as figure 1 for the transition lines. 



investigated which allows us to identify the presence of a mixed phase AF+SG. 
This mixed phase corresponds to a replica symmetry breaking (RSB) SG with 
™P 7^ (p = o.,b) [23 ■ The emergence of each type of solution depends on the 
relationship among parameters g, ( Jo)~^ (the degree of frustration) and F given 
in units of J. 

Therefore, we can build, in the beginning, two kinds of phase diagrams T (T 
is the temperature) versus: (a) g {g is the strength of intrasite pairing interac- 
tion) with Jo and F kept independents; (b) F with g and Jo kept independents. 
The first phase diagram can show directly the competition among AF, SG and 
PAIR phases. The second one gives more precise information about the role of 
the transverse field on the transition lines present in the problem. 

In Fig. 1, we show the results for T versus g for Jo = 1.5J with F = 
and J. Therefore, in the region of small g in Fig. 1, it is quite clear that the 
parameter Jo is related with the presence of magnetic solutions. Firstly, the AF 
solution appears below Tjv. Then, when temperature is decreased, there is the 
onset of a mixed phase AF+SG at Tf. Finally, at lower temperature, there is 
a transition from AF+SG to a SG phase at Tg. The three magnetic transition 
lines mentioned above are second order. For large g, the solutions found for the 
order parameters indicate the existence of the PAIR phase EH EI in which 
there is pairing formation in both sublattices. The role of the transverse field F 
is also clear in this particular phase diagram. The field decreases simultaneously 
the magnetic transition temperatures T/v, Tf and Tg. It also displaces the PAIR 
transition line Ti{g) in the sense that it is necessary to increase the parameter 
g to find a PAIR solution in the problem. 

Fig. 2 shows the previous phase diagrams when the degree of frustration 
is decreased (Jo = 1.7J) with the transverse field kept F = and J. We can 
compare the results in the small and large g region in this figure with the phase 
diagram given in Fig. 1. The conclusion is direct, in the small g region, there are 
competing effects due to Jq and F . The decrease of the degree of frustration 
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0,5 1 1.5 2 2.5 0.5 1 1.5 2 2,5 

r/J r/J 



Figure 3: Phase diagrams as a function of T/J and T / J for Jo — 1.5 J and 
several fixed values of g/J: (a) g = 0, {b) g — 6. 5 J, (c) g = 8.0 J and (d) 
g = 9.0 J. It is used the same convention as figure 1 for the transition lines. 



enhances T^y. However, it decreases Tj and Tg as well. On the other hand, the 
whole set of magnetic transition temperatures Tjv, T/ and Tg decreases with T . 
The superposition of both effects is responsible by the suppression of the SG 
phase in Fig. 2. The line transition Ti{g) is not affected by the change of the 
degree of frustration (Jo)~^ while the transverse field F has the same role as 
before, it displaces Ti{g). 

The numerical analyses indicate that the PAIR transition line Ti{g) is more 
complicated than the magnetic ones. It is a second order phase transition at 
higher temperatures and a first one at lower temperatures, where there are mul- 
tiple PAIR solutions. In this case, the stable solution minimazes the thermody- 
namic potential. This same criterion is used to obtain the first-order boundary. 
The transition lines can be also analysed by performing a Landau expansion 
of the thermodynamic potential in powers of the order parameters {qa, Qb, ^a, 
nib, Va and rjt). We can explore the symmetry of the parameters: q — qa — qt, 
X = Xa = Xb, V — Va — Vb and ma = —mi,. Equation H20|) is expanded in pow- 
ers of q, T] and / = {ma — mb)/2 {I is the antiferromagnetic order parameter), 
while xiQi^iV) is given by the saddle-point equation H25|l . After some lengthy 
calculations, the Landau expansion of the thermodynamic potential is: 

+B2q^ + C2Tf + C^jf (28) 
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with 



A2 = 4/3Jo(l - /3JoXo)P (29) 
2! 



B, = -^+p\j^xl (30) 



G2 = f3g^!^, (31) 



C4 = i7^(— i^oXoX2+3-ifo) (32) 



where ii'o = 1 + cosh a/Aq 

Xo = ^ I D^ r'""Zr\ (33) 




= ~^^^^^°/^° (34) 

1 + PU^XI - / i?^^4(£25^ _ £ml^)/i^o] 

and Aq = 2/3^J^xo£,^ + /3^r^. The tricritical point Ttrict can be obtained from 
Eqs. H31I32|I which show that the transverse field T affects the location of the 
Ttrict- Actually, it moves upwards Ttrict (see Figs. 1 and 2) 20', 24 . 

In Fig. 3, we show the phase diagram T versus T for Jo = 1.5J and g — 
0, 6.5J, 8J and 9J. The case in which there is no pairing coupling is shown 
in Fig. 3-a. In fact, this situation has been studied in Ref. (25, where the 
increase of F leads the transition temperatures T/v, Tf and Tg towards their 
respective QCP's. In particular, the critical transverse field for AF transition 
can be obtained analytically by expanding the sublattice magnetization rup (see 
Ref. [H]). The ordering Tg < Tf < Tm is kept when F increases, which is 
the reason why the transition temperatures Tjv, Tf and Tg are simultaneously 
depressed in Figs. 1 and 2. The increase of g (see Figs, lb, Ic, Id) allows the 
existence of a PAIR solution which depends on, as discussed in Ref. [23 j of 
the ratio V / g as well as the existence of a magnetic solution depends on F/Jq 
in the present work. From this view point, the role of F in Figs. 1 and 2 is 
clearly confirmed in Fig. 3, it tends to suppress any phase which is appearing 
as solution in the problem. Actually, the displacement of the PAIR phase in 
Figs. 1 and 2, when F increases, reflects this effect. Moreover, the increase of g 
also leads (if Jo is kept constant) the PAIR phase to become dominant. It also 
moves upward the tricritical point Ttrict as in Ref. |24| . 

The information contained in Figs. 1-3 can be displayed in a more adequate 
format if we assume a relationship among the parameters F, Jq and g. This 
kind of procedure has already been adopted in Ref. [23 j which is ultimately 
justified by the fact that both RKKY and the pairing interaction in Eq. 
are originated from the same source ^Hl- In the present case, F would have the 
equivalent role of spin flipping part of the Heisenberg model [21]. Besides, it is 
a more convenient format to compare with experimental results. Therefore, we 
assume the following relationship: 

F — aig + Si (35) 
Jo = a2g + 52. (36) 
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Figure 4: Phase diagram T/ J versus g/ J builds for the relation Jq = 3.50—0.715 
and r = 0. The dashed line indicates a first-order phase transition while the 
full lines indicate second-order phase transitions. 



The complicated interplay between Jq and T can be adjusted by the factors 
oil, 012, 5i and 82 in Eqs. (|!l5|l - p(H) . We choose a2 < 0, which means that 
the degree of frustration enhances with the increases of the pairing coupling g. 
The factor 82 is adjusted to guarantee the AF coupling and, with the remaining 
factors, also to maintain the transitions and the tricritical point located at the 
same scale as Figs. 1-3. 

Fig. 4 shows the phase diagram T versus g with F = 0. For that case, 
the behaviour of transition temperatures is obtained from the solution of order 
parameters (Eqs. ^^-^E^) together with the AT hne (Eqs. ^^-^7^) using 
only Eq. H3()|l . For small g (small degree of frustration), there is a transition 
from paramagnetism (NP) to AF phase. Consequently, the Neel temperature 
decreases when g is increased. Then, in a interval of g, there is direct transition 
from NP to SG phase (in this case Tg — Tj). For large g the PAIR phase is 
completely dominant. At temperature T < Tg, the situation is richer when g 
increases as consequence of the RS lack of stability. The solutions found are 
AF at small g and SG at some interval of g as before. However, in a very small 
range of g, a mixed phase AF-I-SG intermediated between AF and the SG phases 
appears. Then, after the sequence of second order transitions AF-AF-I-SG-SG, 
there is a first order boundary between SG and the PAIR phase. The location 
of the tricritical point Ttrict is not changed as compared with phase diagrams 
shown in Figs. 1 and 2 when F = 0. 

Figs. 5-6 show the solution for the order parameters when the transverse 
field r is also tunned by Eq. (jSHJ- The numerical factors ai and 61 can be 
used to adjust the strength of the transverse field as g increases. In Fig. 5, 
r affects the transition line Tn and particularly Tg, which moves downwards 
when g increases. The sequence of phases at lower temperature is preserved 
as in Fig. 4. However, F is not strong enough to lead Tg to a QCP. On the 
other hand, the behaviour of the PAIR phase boundary Ti{g) is affected as 
in Figs. 1-3 in which Ttrict moves upward when F increases. In particular, it 
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Figure 5: Phase diagram T/J versus g/J obtained from a relationship among 
Jo/ J, T/J and g/J (Jo = 2.77 - 0.24^ and V = 1.30 + O.O95). It is used the 
same convention as figure 1 for the transition hnes. 



is possible to find one metastable SG solution into the PAIR phase below Tg 
which keeps going downward. Fig. 6 displays the situation where ai and 5i are 
adjusted to enhance the strength of F as compared with Fig 5. For that case, 
the spin flipping induced by F is strong enough to lead Tg to a QCP while the 
Ttrict is obtained at a larger value than before. For both cases of Figs 5 and 6, 
the tunning of F still preserves the sequence of phases AF-AF+SG-SG at low 
temperature likewise the case F = for a certain range of g. 

4 Conclusions 

In the present paper, we have analysed the competition among antiferromag- 
netism ( AF) , spin glass (SG) and pairing formation phase (PAIR) in the presence 
of a quantum tunneling mechanism. The two-sublattice model used is composed 
by a Gaussian random interlattice Ising interaction (with mean Jo and standard 
deviation J) j26| . an intralattice pairing interaction with an applied transverse 
field F. The partition function is calculated in the functional integral formalism 
in which the spin operators are given by bilinear combinations of Grassmann 
fields |19l 1201 I24j . The saddle-point Grand Canonical potential is obtained 
within of static approximation (SA), the replica symmetry (RS) ansatz and in 
the half-filling. Particularly, the use of SA is justified because our main interest 
is to study in detail the phase boundaries among AF, SG and PAIR phases 
when the spin fiipping is activated by a transverse field F. 

In the mean field theory presented, the phase transitions of the fermionic 
system defined in Eq. appear related with pairing and magnetic internal 
fields for each sublattice p. The magnetic one hp {p — A, B) has a random and 
AF components. In particular, the AF part of hp depends on the sublattice 
magnetization mpi as well as the random part is associated with the replica 
non-diagonal SG order parameter qpi and also with Xpi — Qpi — Qp' {qpi is the 
replica diagonal SG order parameter), where p p' . In contrast, the pairing 
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Figure 6: Phase diagram T/J versus g/J for the relations Jq = 2.00 — 0.322(7 
and r — 1.85 + 0.09(7. The same convention as figure 1 is used for the transition 
Unes. 



internal field applied in the sublattice p depends on the PAIR order parameter 
rjp. Furthermore, there is the presence of F which tunes the spin flipping and, 
hence tends to suppress any kind of magnetic phase. The pairing formation is 
also affected, as it can be clearly seen in Fig. 3. 

The solutions for qp, x^, rup and rjp (PAIR order parameter) are located in a 
parameter space given (in units of J ) by Jqi T and g (the strength of the pairing 
interaction). Figures 1, 2 and 3 show the phase diagrams for several cuts in the 
previous space. Thus, it is possible to identify how each parameter can favour 
one particular solution for a given temperature. Another important point is to 
locate the Almeida- Thouless line (T/) in such space. To take a typical case, in 
Fig. 2, we present the phase diagram temperature versus g for a Jq = 1.7 J and 
F = and J. For lower g, the magnetic solutions are dominant with a sequence 
of second order phase transitions AF, a mixed phase AF+SG and SG at lower 
temperature. For larger g, the local pairing is dominant. Actually, the PAIR 
phase boundary Ti{g) has a complex nature with the presence of a tricritical 
point Ttrict which is quite dependent on the transverse field F (see Eqs. H31I32II ). 

In the phase diagrams 4, 5 and 6, we propose a relationship between the 
parameters Jq and F with g (see Eqs. (|35l36|l l based on the original derivation 
of the model given in Eq. iQJ (see Ref. ) • This procedure allows to compare 
our results with the phase boundaries found in experimental phase diagrams. 
For instance, there are some similarities between the experimental situation for 
Yi_a; CaxBa2CuzO(i and the phase diagram shown in Fig. 5 as well as between 
Ui^xLaxPd2Al^ and the one shown in Fig. 6, if it is possible to associate 
the doping in those physical systems with the parameter g. For the first case, 
similarities such as the sequence of phases and, in particular, the presence of SG 
mixed with an AF background. In terms of the present model, the decrease of 
experimental SG temperature transition Tg could be explained by the presence 
of quantum spin flipping mechanism which is not strong to lead Tg towards a 
QCP. For the second one, there are also similarities not only between the phase 
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boundaries, but also with the behaviour of Tg which is depressed to a QCP. 

To conclude, in this work we studied the thermodynamics of the model given 
in Eq. Our goal is to obtain the corresponding phase boundaries and, then 
to mimic the global phase diagram of physical systems as Yi^xCaxBa2Cu^OQ 
and Ui-xLaxPd2Al^. As last remarks it should be noticed the role of Xp, which 
carries the effects of disorder even at T > Tg, to determine the PAIR phase 
boundary Ti{g). In the case T > Tg, Xp is equal to the replica diagonal SG 
order parameter g^. In Fig. 6, because of the presence of QCP, this identity is 
true for the entire range of temperature after the QCP. It is well known that 
Qp can be written in terms of the site occupation in the sublattices At the 
same time, the nature of Ti{g) depends deeply on the transverse field F. This 
arises the question how the phase boundaries in the present problem would be 
affected by the interplay between the chemical potential n and F in situations 
like those shown in Figs. 5 and 6. It is also important to remark that the precise 
location of the phase boundaries below Tf needs RSB spin glass solutions, as 
for instance the boundary between SG and mixed phase and the first-order 
transition between SG (or mixed phase) and PAIR phase |23- The study of the 
situation where 7^ and the implementation of RSB will be object of a future 
work. 
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